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Ice-Water  Partition  Coefficients  for  RDX  and  TNT 

SUS,\N  TAYLOR 


INTRODUCTION 

The  purpose  of  this  study  was  to  measure  the 
partition  coefficients  (between  water  and  ice)  of 
RDX  and  TNT  and  to  determine  how  rapidly  these 
substances  can  diffuse  away  from  a  freezing  front. 
Partition  coefficient  is  used  here  in  a  nonrigorous 
sense  and  describes  how  much  TNT  or  RDX  is  in 
the  liquid  versus  the  ice  phase  at  nonequilibrium 
conditions.  Studies  of  this  kind  have  been  done 
using  brine  (Cox  and  Weeks  1975),  but  no  data  exist 
on  the  partition  coefficien*-  between  ice  and  water 
for  organic  molecules. 

This  work  was  conducted  as  part  of  a  research 
program  investigating  the  use  of  artificial  freezing 
for  decontamination  of  soils.  Since  RDX  and  TNT 
are  being  used  to  see  if  freezing  can  move  organics 
in  saturated  soils,  knowledge  of  their  behavior,  i.e., 
diffusivity  and  partition  coefficient,  would  help 
determine  the  conditions  under  which  they,  and 
the  volatile  organics  of  interest  for  hazardous  waste 
cleanup,  might  be  moved.  An  ice-water  partition 
experiment  using  RDX  and  TNT  was  conducted  to 
determine  the  efficiency  writh  which  the  formation 
of  ice  excludes  nonvolatile  organic  compounds.  An 
upper  limit  for  the  freezing  rate  needed  to  exclude 
large,  nonvolatile  organic  molecules,  i.e.,  TNT,  RDX 
etc.,  from  the  ice  structure  was  calculated  for  this 
experiment. 

EXPERIMENTAL  METHODS 

A  1-mg/L  solution  of  RDX  and  TNT  in  reagent 
grade  water  was  prepared  using  production  grade 
RDX  and  TNT  purchased  from  Eastman  Kodak. 
The  solution  was  made  up  in  two  20-L  jugs  and  was 
mixed  continuously  for  two  days.  The  jugs  were 
covered  with  tinfoil  to  exclude  light.  The  solution 
was  then  transferred  to  a  30-cm-wide  by  40-cm- 
long  by  30-cm-high  reservoir  and  subjected  to  a 


0.367°C/cm  temperature  gradient  (ll°C/30  cm). 
The  base  of  the  Plexiglas  container  was  maintained 
at  l^C  and  the  freezing  plate,  which  also  served  as 
the  top  of  the  container  and  was  in  contact  with  the 
solution,  was  kept  at  -10°C.  Two  holes  at  the  base 
of  the  container  served  as  overflow  condu  i  ts  for  so¬ 
lution  pushed  out  during  ice  formation.  The  sides 
of  the  container  were  insulated  with  polystyrene 
insulation  board.  This  experimental  setup  (see  Fig. 
1)  resulted  in  a  horizontal  ice/ water  boundary,  in¬ 
dicative  of  unidirectional  freezing,  in  which  the  ice 
grew  from  the  top  downward. 

The  freezing  rate  was  not  constant,  being  faster 
td  the  beginning  than  at  the  end  of  the  experiment. 
The  ice  thickness  was  periodically  marked  on  the 
Plexiglas  container  and  an  estimate  of  the  freezing 
rate  was  obtained  by  plotting  the  ice  thickness 
versus  time.  Unfortimately,  the  ice  growth  was  not 
monitored  until  the  second  day  (when  I  discovered 
that  the  data  logger  connected  to  the  thermocouples 
was  not  working  properly),  so  information  on  the 
fastest  ice  growth  is  not  available.  Based  upon 
information  derived  from  ice  growth  in  brines,*  a 
3rd  order  polynomial  cvuve  was  fit  through  the 
points  and  the  origin  (Fig.  2)  to  project  ice  growth 
on  the  first  day. 


40  cm 


Figure  1 .  Experimental  setup. 


Tersonal  communication  with  Dr.  Gordon  Cox,  CRREL,  June 
1988. 


Figure  2.  Plot  of  ice  thickness  versus  time. 

Roughly  one- third  (10.1  cm)  of  the  reservoir  was 
frozen  by  the  end  of  the  experiment.  A  3-cm  section 
was  removed  from  each  edge  of  the  ice  block  to 
minimize  any  effects  the  side  wall  might  have  had 
on  the  growth  and  hence  the  concentration  of  the 
ice.  The  central  section  was  used  for  the  analyses. 
Figure  3  illustrates  the  manner  in  which  the  origi 
nal  ice  block  was  subdivided.  To  obtain  a  general 
idea  of  the  distribution  of  RDX  and  TNT  in  the  ice, 
1  -cm-thick  sections  of  block  H  were  cut  with  a  band 


!  Table  1.  Concentrations  measured 

for  known  standards  and  explo- 
I  sive  controls.  The  variation  between 

d  the  known  value  and  those  measured 

— i  is  about  5%.  The  analyses  of  the  ini- 

-j  tial  and  final  solutions  are  also  tabu¬ 

lated. 


RDX 

TNT 

<\iglU 

Explosive  standard 

2006 

365 

Explosiv'e  standard 

1915 

337 

(measured) 

2095 

354 

2028 

342 

Explosive  control 

220 

207 

Explosive  control 

236 

193 

(measured) 

227 

186 

Initial  solution 

1181 

1160 

(measured) 

1223 

1192 

Final  solution 

1831 

1696 

(measured) 

1831 

1670 

1831 

1683 

00 

1676 

and  a  Cole  Parmer  chart  recorder.  An  LC-18  col¬ 


saw  and  analyzed.  After  the  analysis  of  these  initial 
samples  indicated  that  the  results  looked  promis¬ 
ing,  a  more  careful  sampling  of  the  ice  block  was 
obtained  by  shaving  blocks  B  and  E  with  a  mi¬ 
crotome.  A  series  of  very  thin  layers  was  removed 
from  the  ice  and  placed  in  bottles  of  known  weight. 
Weights  of  the  sample  were  obtained  and  just 
enough  sample  was  coUected  to  fill  the  2-mL  HPLC 
autosampler  vials.  Samples  of  the  starting  solution 
and  the  final  solution  were  also  analyzed. 

All  the  chemical  analyses  were  conducted  using 
a  Reversed-Phase,  High  Performance  Liquid  Chro¬ 
matography  (RP-HPLC)  instrument  composed  of  a 
Spectra  Physics  detector  (SP8490)  and  pump 
(SP881 0)  coupled  to  a  Dynatech  autosampler  (model 
LC  241),  a  Hewlett-Pack2u-d  integrator  (HP3393A) 


umn  (Supelco,  Inc.)  was  eluted  with  1:1 
methanol-water  at  a  flow  rate  of  1.5  mL/ min.  The 
variable  wavelength  detector  was  set  at  254  nm. 
Retention  times  for  RDX  and  TNT  were  4.3  aiid  9.4 
minutes  respectively.  Peak  heights  of  2  cm  for  RDX 
and  5  mm  for  TNT  at  a  range  setting  of  0.CX)1  ab¬ 
sorbance  units,  full  scale,  corresponded  to  the  de¬ 
tection  limits  and  could  be  measured  accurately  by 
hand.  At  the  low  concentrations  of  interest  to  this 
study,  manual  peak  height  measurements  proved 
to  be  more  reliable  than  automated  peak  areas  from 
the  integrator. 

The  melt  water  from  the  ice,  the  starting  solution 
and  the  unfrozen  final  solution  contained  no  solids 
and  were  run  on  the  HPLC  without  diluting  or 
filtering.  Table  1  lists  the  concentrations  measured 
for  the  starting  solution  (average  1.202  mg/L  for 
RDX  and  1.176  mg/L  for  TNT)  and  for  the  final 


Ice/Wofer  Interface 


lO.l  err 


3  cm  Removed 
(oil  around) 


40  cm 


Ice/Plote  Interface 


Figure  3.  Sectioning  of  the  ice  block. 


solution  (average  1.829  mg/L  for  RDX  and  1.681 
mg/L  for  TNT)  using  the  HPLC.  Explosive  stan¬ 
dards  and  an  explosive  control  were  analyzed  along 
with  the  samples  to  ensure  that  the  iiwtrument  was 
operating  properly. 

As  the  solubilities  of  RDX  and  TNT  in  water  are 
low — ^30  mg/L  for  TNT,  40  mg/L  for  RDX  (Leggett 
1985) — an  origina]  concentration  of  1  mg/L  was 
chosen  to  ensure  that  any  increase  in  TNT  and  RDX 


I 


concentration  in  the  water,  resulting 
from  their  exclusion  from  the  ice,  would 
not  result  in  precipitation  of  the  solutes. 
The  concentration  selected  also  made 
sample  dilution  unneccessary  when 
RDX  and  TNT  were  determined  by  RP- 
HPLC. 


RESULTS  AND  DISCUSSION 

The  majority  of  the  ice  block  was  co¬ 
lumnar  growth  icc,  with  a  7-mm-thick 
transition  zone  at  the  freezing  plate 
interface  (Fig.  4).  The  columnar  ice  had 
fairly  equidimensional  grains  (Fig.  5), 
whereas  the  transition  zone  had  grains 
that  had  unequal  dimensions  (Fig.  6). 
Many  small  bubbles  within  the  transi¬ 
tion  ice  made  this  region  opaque.  Tubu¬ 
lar,  vertical  bubbles  were  found  in  the 
columnar  ice  (see  Fig.  7).  The  different 
bubble  patterns  are  a  function  of  the 
rates  of  freezing  and  the  rates  of  diffu¬ 
sion  of  gas  in  water  (Halde  1979).  At 
high  freezing  rates,  numerous  small 
bubbles  form,  whereas  lower  freezing 
rates  give  the  dissolved  gases  time  to 
diffuse  to  existing  bubbles.  At  very  slow 
freezing  rates,  the  dissolved  gas  has 


figure  4.  Vertical  thin  section  of  the  ice  block. 


Figures.  Horizontal  section  through  the  columnar  ice. 
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time  to  diffuse  a  wav  from  the  freezing  front  and  no 
bubbles  form.  This  was  the  case  for  the  bottom 
mm  of  the  ice. 

Preliminary  chemical  analyses  using  1-cm-thick 
sections  (120-mL  volumes)  indicated  that  only  llie 
ice  at  the  two  interfaces — the  ice/ freezing  plate  in¬ 
terface  and  the  ice/water  interface — had  any 
measurable  RDX  and  TNT.  Note  that  the  method 


detection  limit  is  12  |Jg/L  for 
RDX  and  2.6  |ag/L  for  TNT 
(Jenkins  et  al.  1988). 

A  detailed  chemical  analy¬ 
sis  of  the  ice  at  the  two  inter¬ 
faces  was  then  made.  Small 
samples  (2  mL)  v\'ereobtained 
by  microtome  sectioning.  Fig¬ 
ure  8  is  a  plot  of  RDX  and 
TNT  concentration  versus  po- 
sition  of  the  ice  below  the 
freezing  plate  obtained  from 
the  latter  analysis.  Concen¬ 
trations  were  below  the  de¬ 
tection  limit  adjacent  to  the 
freezing  plate,  increased  to  a 
maximum  at  about  4  mm,  fell 
below'  the  detection  limits  at 
1  cm  and  then  increased  right 
at  the  ice/ water  interface. 

The  low  concentration  seen 
in  the  first  millimeter  is  inter¬ 
preted  as  resulting  from  a  low 
freezing  rate  experienced  briefly  between  the  time 
the  freezing  pla  te  was  turned  on  and  when  it  equili¬ 
brated  at  -1()°C.  Once  the  plate  w'as  cold,  however, 
the  freezing  rate  increased,  entrapping  up  to  15% 
TNT  and  RDX.  The  progressive  decrease  in  concen¬ 
tration  after  the  maximum  reflects  a  gradual  de¬ 
crease  in  the  freezing  rate  brought  about  by  the  ice 
insulating  the  freezing  plate. 


Figure  7.  Appearance  of  the  ice  during  the  course  of  the  experiment. 


=  changing  concentrarion 
of  the  water 

C  =  initial  water  concentra- 

C1 

tion 

=  average  concentration 
of  the  ice 
h  -  tank  depth 
t  =  ice  thickness. 


Coj^cen^rason  ug/L) 

Figure  S.  Meustiml  concentratums  ofRDX  and  TNT  as 
depth  loithin  the  ice  block. 

The  measured  concentration  (60  g/L  for  each 
compound)  in  the  ice  at  the  ice/water  interface  is 
attributed  to  the  adherence  of  a  small  amount  of  the 
final  liquid  to  the  ice  surface.  This  interpretation  is 
supported  by  my  having  observed  drops  of  water 
on  the  surface  of  the  ice  as  the  water  was  drained 
from  the  base  uf  the  tank.  As  the  ice  was  kept  below 
freezing,  the  drops  were  not  produced  by  melting. 
Originally,  I  thought  the  spike  in  concentration  at 
the  ice/water  interface  was  attributable  to  an  ad¬ 
hered  monolayer.  The  120  g/L  (60  g/L  of  each 
compound)  observed  is,  however,  less  than  the 
2.35  mg/L  piedicted  for  a  monolayer  (see  Appen¬ 
dix  A). 

In  order  to  calculate  the  effective  partition  coef¬ 
ficient,  C./ C^^,,  not  only  is  the  concentration  of 
TrjT  and  RDX  in  the  ice  required  but  that  in  the 
remaining  water  is  also  needed.  Since  direct  sam¬ 
pling  of  water  in  the  reservoir  was  not  possible 
during  the  course  of  the  experiment,  the  initial  and 
final  concentrations  of  TNT  and  RDX  in  the  water 


When  the  concentration  of  TNT 
and  RDX  in  the  ice  was  below  the  detec¬ 
tion  limit,  the  detection  limit  for  each 
organic  was  used  as  the  contaminant 
concentration  in  the  ice.  The  calculated 
final  concentration  for  RDX  compared 
well  with  the  value  measured  with  the 
RP-HPLC(1820  g/L  calculated,  1829 
g/Lmeasured)butforTNT  there  was 
r .  ■  -  a  6%  difference  between  the  calculated 

; _ '  and  the  measured  value  (1787  g/L 

ooc  calculated,  1681  g/L  measured).  Be¬ 
cause  the  concentration  of  TNT  and 
hinction  of  so  much  larger  in  the  water  than 

the  ice,  the  calculated,  final  concentra¬ 
tion  is  not  sensitive  to  the  value  chosen 
for  the  below-detection  ice  concentra¬ 
tion.  In  the  future  the  liquid  expelled  out  the  tube 
will  be  sampled  with  time  and  analyzed. 

From  the  C,  measured  and  the  C  calculated,  the 
effective  partition  coefficient  was  determined 
(Table  2).  Following  Cox  and  Weeks  (1975),  the 
equation  of  Burton  et  al.  (1953)  was  used  to  relate 
the  effective  partition  coefficient  (X^^^)  to  the  "theo¬ 
retical"  partition  coefficient  experienced  at  an  ice 
growth  velocity  v  =  0 

Y  -  (j) 

K  pff - —  -  — 

tC  j-  ti 

+  n-^o)e 

where  X^jj  =  effective  partition  coefficient 

X  =  partition  coefficient  v.’hen  freez- 
ing  front  is  static,  v  =  0 
V  =  velocity  of  freezing  front 
D  =  diffusion  coefficient  of  RDX  and 
TNT  in  water 


and  the  measured  TNT  and  RDX  concentrations  §  =  thickness  of  the  boundary  layer, 

for  each  sampled  ice  layer  were  used  to  calculate 

the  changing  TNT  and  RDX  concentrations  in  the  By  rearranging  eq  2  in  linear  form 
bulk  water.  The  following  equation  from  Cox  and 
Weeks  (1975)  was  used. 

1 

C  =lCh-C.t]/h-t  (1) 

W  O  I 
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Tabk2,  Calculated  values  of  for  TNT  and  Table  3.  Organic  substances  whose  diffu- 

RDX  at  different  ice  layers.  The  ice  growth  veloc-  sivities  and  molecular  weights  were  used 

ity  and  the  position  of  each  layer  within  the  ice  in  Figure  10. 

block  are  also  tabulated. 


Mt>hviilar*' 

/a- 

Oi/fiisii'i(ics* 

u’l'iy/it 

t'clocity 

Sulvlaiicc 

(D  XI0.S) 
a 

(a.iii.ii.) 

VM 

(dll' 

(llll/s) 

TNTK. 

RDX 

4Tr 

lYr 

Ljctow 

0.43 

342.3 

18..30 

!).) 

1.08x10-* 

2.54x10'* 

4.57x10'* 

Sucrose 

0.45 

342.3 

18.50 

0.2 

l.lMxlO-* 

1.52x10'- 

1.65x10“" 

Caffeine 

0.63 

194.2 

13.94 

0.3 

1.08x10-* 

7.41x10'- 

6.60xUr" 

Mannitol 

0.58 

182.18 

13.50 

0.4 

1.04xU)-* 

1.16x10“' 

1.07x10“' 

Glucose 

0.6 

180.16 

13.42 

0.3 

l.lb.xlO-* 

1.00x10“' 

9.82x10“" 

Tyrogallol 

0.7 

126.11 

11.23 

O.h 

1.08x10-* 

9.01x10'^ 

8.57x10'" 

Rosorcinol 

O.S 

110.11 

10.49 

0.7 

1.01x10-* 

6.65x10'- 

6.51x10'" 

Hydroquinone 

0.77 

IIO.I! 

10.49 

0.8 

1.02x10-* 

3.81x1  O'" 

4.05x10'" 

Phenol 

0.84 

94.11 

9.70 

0.0 

1. 00x10“* 

1.53x10'" 

1.62x10'" 

Pyridine 

0.76 

79.1 

8.89 

1 

1.02x10“* 

5.43x10'* 

9.66x10'* 

.  - 

■) 

0.09x10-* 

1.59x10'* 

9.32x10'* 

*T}iibcKleaiiv(1979,  p.  462-463) 

3 

7.70x10-* 

1.53x10* 

8.99x10'* 

*'■  Atomic  mass  units. 

4 

7.10x10-* 

1. 47x1  O'* 

8.65x10'* 

5 

b.40xl0-* 

1.41x10“* 

8.31x10'* 

h 

4.74x10  * 

1. 36x10“* 

7.79x10'* 

7 

3.79x10-* 

1.30x10'* 

7.64x10'* 

8 

3,13xUr* 

1.24x10'* 

7.30x10'* 

2.63x10“* 

1.18x10'* 

6.96x10'* 

U) 

2.86x10“* 

1.12x10'* 

6.62x10'* 

one  can  see  that  by  plotting  the  natural  log  of  K  limits,  the  values  obtained  for  being  very  sensi- 

against  the  ice  growth  velocity,  v,  the  slope  of  the  tive  to  the  value  chosen  for  Cj  when  the  concentra- 

resulting  line  is  the  value  of  S/  D  with  being  the  lion  of  the  organics  in  the  ice  was  below  the  detec- 

intercept  (Cox  and  Weeks  1975).  Figure  9  is  such  a  tion  limit. 

plot  and  yields  an  upper  limit  for  K^,  the  partition  Since  diffusivities  of  molecules  in  water  depend 

coefficient  at  u  =  0,  of  6.06  x  10'^  for  W)X  and  upon  the  size  of  the  molecule,  a  plot  of  measured 

1.07  X  ICT^  for  TNT.  These  values  are  only  upper  diffusivities  of  various  organic  molecules  (Table  3) 


1  ’  '  r--|— 
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Fiyjiiv  10.  Mcnf  lin'd  diffiifiivifii'^  often  on^iinicf 
vf  rhaquiiiv  root  of  thciniioiccular 'aright.  Sluuicd 
atvti  ivpn':<oiit<  one  standard  deviatton  from  the 
best  tit  line. 


the  square  root  of  their  molecular  weights 
gi\'es  a  linear  relationship’*'  (Fig.  10).  Plotting  TN’T 
and  RDX  on  this  graph  gives  an  appro.ximatediffu' 
Sion  coefficient  for  RDX  of  5.8  X  10'^’  cm  /s  and  for 
TNT  of  t.7  X  lO’*’  cm  /s  .  By  multiplying  the  values 
obtained  for  5/D  from  Figure  9  (3736  s/cm  for 
RDX,  4396  s /cm  for  TNT)  by  their  respective  values 
for  D,  5  was  estimated  to  be  2.2  x  10'-  cm  for  RDX 
and  2.5  x  10"-  cm  for  TNT. 


CONCLUSION  AND 
FUTURE  RESEARCH 

In  conclusion,  I  think  that  TNT  and  RDX  are 
suitable  for  modeling  the  movement  of  organics 
upon  freezing,  as  they  are  excluded  from  the  ice 
and  are  not  volatile  \n  upper  limit  for  the  partition 
coefficient,  K  was  calculated  to  be  6.06  x  10'^  for 
RDX  and  1 .07  x  10"^  for  TNT.  These  values  were 
obtained  by  calculating  the  effective  partition  coef¬ 
ficients  (concentration  of  TNT  and  RDX  in  frozen 
layers  versus  the  changing  liquid  composition) 
and  plotting  these  values  against  the  ice  growth 
rate.  The  thickness  of  the  boundary  layer  was  cal¬ 
culated  to  be  about  2.5  x  10"^  cm. 


^Personal  communication  vvitli  Daniel  Leggett,  CRREL,  June 
1988. 


Freezing  excludes  TNT  and  RDX  very  efficiently 
at  freezing  rates  up  to  9  x  10’'^cm/s  (3.1  in. /day);  at 
higher  freezing  rates  up  to  15%  of  the  solute  was 
entrapped.  TNT  and  RDX  behaved  similarly  in  this 
experimen  t,  which  is  not  surprising  given  the  close¬ 
ness  of  their  molecular  weights  (227  for  TNT,  222 
for  RDX).  These  compounds  should  be  shielded 
from  light  to  avoid  photodegradation.  The  use  of  a 
microtome  as  a  sampling  tool  was  found  to  provide 
a  detail'd  record  of  the  distribution  of  these  mole¬ 
cules  as  a  function  of  freezing  rate. 

In  future  experiments,  soil  or  a  particulate 
medium  (glass  beads)  w'ill  be  used  to  study  the 
effect  particles  have  on  the  concentration  and  dis¬ 
tribution  of  nonvolatile  organics.  In  my  opinion, 
the  variables  of  greatest  importance  for  determin¬ 
ing  the  distribution  of  the  organics  are:  tortuosity 
(path  length),  chemical  activity  of  the  particulate 
material  (both  organic  and  inorganic)  and  the  freez¬ 
ing  rate.  Glass  spheres  or  coarse  quartz  sand  should 
give  information  on  the  effect  of  tortuosity  on  the 
movement  of  organics.  The  added  complication  of 
surface  chemistry  would  be  studied  using  soils. 
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APPENDIX  A:  CALCULATION  OF  CONCENTRATION  OF  MONOLAYER 


The  concentration  of  the  monolayer  was  calcu¬ 
lated  by  dividing  the  area  of  the  ice  surface  (80.1 
cm^)  by  the  surface  area  of  each  compound.  From 
the  number  of  TNT  or  RDX  molecules,  the  number 
of  moles  and  then  the  number  of  grams  was  calcu¬ 
lated.  The  gram  weight  of  the  monolayer  was  then 
divided  by  the  volume  of  the  sample.  The  surface 
areas  of  TT^  and  RDX  were  estimated  using  the 
following  equation  from  Bansal  et  al.  (1982)  (Leggett 
1985): 


A  =  1.091x10-’MMx102VN]^''^  (Al) 

where  A  =  area  in  square  centimeters 
M  =  molecular  weight 
N  =  Avogadro's  number  =  6.023  x  10^ 
molecules/mole. 

Equation  Al  yields  a  surface  area  for  the  TNT 
molecule  of  5.7  x  lO'^^cm^  (5.6  x  10"'^cm^  for  RDX). 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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